Quantum damped oscillator I: 
dissipation and resonances 



in 
o 
o 

(N 



c3 



Dariusz Chrusciriski and Jacek Jurkowski 
Institute of Physics, Nicolaus Copernicus University 
ul. Grudziadzka 5/7, 87-100 Torufi, Poland 



Abstract 



Quantization of a damped harmonic oscillator leads to so called Bateman's dual 
system. The corresponding Bateman's Hamiltonian, being a self-adjoint operator, 
displays the discrete family of complex eigenvalues. We show that they correspond 
to the poles of energy eigenvectors and the corresponding resolvent operator when 
continued to the complex energy plane. Therefore, the corresponding generalized 
' eigenvectors may be interpreted as resonant states which are responsible for the 

f*"*". , irreversible quantum dynamics of a damped harmonic oscillator. 

o 
o , 

^ ■ 1 Introduction 

IT) 

The damped harmonic oscillator is one of the simplest quantum systems displaying the 
dissipation of energy. Moreover, it is of great physical importance and has found many 
applications especially in quantum optics. For example it plays a central role in the 
quantum theory of lasers and masers [1, 2, 3]. 

As is well known there is no room for the dissipative phenomena in the stan- 
dard Hilbert space formulation of Quantum Mechanics. The Schrodinger equation 
defines one-parameter unitary group and hence the quantum dynamics is perfectly 
time-reversible. The usual approach to include dissipation is the quantum theory of 
open systems [4, 5, 6, 7]. In this approach the dynamics of a quantum system is no 
longer unitary but it is defined by a semigroup of completely positive maps in the space 
of density operators [8] (for recent reviews see e.g. [9, 10]). 

There is, however, another way to describe dissipative quantum systems based on 
the old idea of Bateman [11]. Bateman has shown that to apply the standard canon- 
ical formalism of classical mechanics to dissipative and non-Hamiltonian systems, one 
can double the numbers of degrees of freedom, so as to deal with an effective iso- 
lated classical Hamiltonian system. The new degrees of freedom may be assumed to 
represent a reservoir. Applying this idea to damped harmonic oscillator one obtains 
a pair of damped oscillators (so called Bateman's dual system): a primary one and 
its time reversed image. The Bateman dual Hamiltonian has been rediscovered by 
Morse and Feshbach [12] and Bopp [13] and the detailed quantum mechanical analysis 
was performed by Feshbach and Tikochinski [14]. The quantum Bateman system was 
then analyzed by many authors (see the detailed historical review [15] with almost 600 
references!). 
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Surprisingly, this system is still worth to study and it shows its new interesting 
features. Recently it was analyzed in [16] in connection with quantum field theory and 
quantum groups (see also [17, 18]). Different approach based on the Chern-Simons 
theory was applied in [19]. In a recent paper [20] a damped oscillator was quantized by 
using Feynman path integral formulation (see also [21]). Moreover, the corresponding 
geometric phase was calculated and found to be directly related to the ground-state en- 
ergy of the standard one-dimensional linear harmonic oscillator. Bateman's system has 
been also studied as a toy model for the recent proposal by 't Hooft about deterministic 
quantum mechanics [22, 23]. 

In the present paper we propose a slightly different approach to the Bateman system. 
The unusual feature of the Bateman Hamiltonian is that being a self-adjoint operator 
it displays a family of complex eigenvalues. We show that these eigenvalues correspond 
to the poles of energy eigenvectors and the corresponding resolvent operator when 
continued to the complex energy plane. The similar analysis for the toy model of a 
quantum damped system was performed in [24, 25]. Eigenvectors corresponding to 
the poles of the resolvent are well known in the scattering theory as resonant states 
[29, 30]. It shows that the appearance of resonances is responsible for the dissipation 
in the Bateman system. Obviously, the time evolution is perfectly reversible when 
considered on the corresponding system Hilbert space H = L 2 (M. 2 ). It is given by the 
1-parameter group of unitary transformations U(t) = e~ lHt . It turns out that there 
are two natural subspaces S± £ H such that U(t) restricted to S± defines only two 
semigroups: U(t > 0) on 5_, and U(t < 0) on <S+. These two semigroups are related 
by the time reversal operator T (see Section 6). 

Our analysis is based on a new representation of the Bateman Hamiltonian, cf. 
Section 4. This representation is directly related to the old observation of Pontriagin 
[31] (see Section 3 for review) that any non-Hamiltonian system of the form 

Xk = X k (xi, . . . ,x N ) , k = 1,2, . . . ,N , (1.1) 

may be treated as a Hamiltonian one in the extended phase-space (xi , . . . ,xn,Pi, ■ ■ ■ ,Pn) 
with the Hamiltonian 

N 

H(x 1 ,...,x N ,p 1 ,...,p N ) = y^p k X k {xi, . . . ,x N ) . (1.2) 

k=l 

Note, that the above Hamiltonian has exactly the form considered by 't Hooft [22]. 

From the mathematical point of view the natural language to analyze the Bateman 
system is the so called rigged Hilbert space approach to quantum mechanics [26, 27, 28]. 
There are two natural rigged Hilbert spaces, or Gel'fand triplets, corresponding to 
subspaces S±. We shall comment on that in Section 8. 

2 Bateman Hamiltonian 

The classical equation of motion for one-dimensional damped oscillator with unit mass 
reads 

x + 27X + kx = , (2.1) 
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where 7 > denotes the damping constant. Introducing Bateman's dual system 

y -2-fy + ny = , 
one may derive booth equations from the following Lagrangian 

L{x, x, y, y) = xy - nxy + -/(xy - xy) . 
Introducing canonical momenta 

Px = y-7y, p y = x + -fx, 

one easily finds the corresponding Hamiltonian 

H(x,y,p x ,p y ) = p x py - j(xp x -yp y ) + u?xy , 

where 
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= \/k - 7 2 



(2.2) 

(2.3) 

(2.4) 

(2.5) 
(2.6) 



Throughout the paper we shall consider the underdamped case, i.e. k > j 2 . 

Now, assuming symmetric Weyl ordering the canonical quantization is straightfor- 
ward and leads to the following self-adjoint operator in the Hilbert space L 2 (IR 2 , dxdy): 



H = Hq + Hi , 



where 
and 

Note, that 



Ho = PxPy + u 2 xy , 
Hi = ~((xp x + p x x) - (ypy + p y y^j 

[H ,H I } = 0. 



(2.7) 

(2.8) 
(2.9) 
(2.10) 



Following Feshbach and Tichochinsky [14] one introduces annihilation and creation 
operators 



A = 



B 



1 



{Px +P y ) ~ iu(x + y) 
{px ~ P y ) ~ iu(x - y) 



They satisfy the standard CCRs 



[i,it] = [B,Bt] = l, 



(2.11) 
(2.12) 

(2.13) 



and all other commutators vanish. It turns out that the transformed Hamiltonian is 
given by (2.7) with 



(2.14) 
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It is easy to see [14, 16] that the dynamical symmetry associated with the Bateman's 
Hamiltonian is that of SU(1, 1). Indeed, constructing the following generators: 

Ji = ^(A*& + AB), (2.15) 

J 2 = % -{ft&-AB), (2.16) 

J 3 = \{tiA + B&), (2.17) 

one easily shows that they satisfy su(l, 1) commutation relations: 

[J1J2] = iJs , [4 U = iJi , [Ji, Js] = iJ2 ■ (2.18) 
Moreover, the following operator 

J = ±(AlA-&B), (2.19) 
defines the corresponding su(l, 1) Casimir operator. One easily shows that 

Jo=\ + J i ~ Ji ~ Ji ■ (2- 2 °) 

It is therefore clear that the Hamiltonian (2.14) can be rewritten in terms of su(l, 1) 
generators as 

Hq = 2HujJ , Hi = 2h-fJ 2 . (2.21) 

The algebraic structure arising in this approach enables one to solve the corresponding 
eigenvalue problem. Let us define two mode eigenvectors \n A ,n B ): 1 

A'A\n A ,n B ) = n A \n A ,n B ) , -B 1 B\n A ,n B ) = n B \n A ,n B ) ■ (2.22) 

It is convenient to introduce 

j = ^( n A-n B ), m=^(n A + n B ), (2.23) 

and to label the corresponding eigenvectors of Jo and J3 by \ j, m) rather than \n A , n B ): 

MiM = j\j,m) , (2.24) 
J 3 \j,m) = (m+^J\j,m). (2.25) 

Clearly, 

j = 0,±± ±1,±|,..., m = \j\,\j\ + l,\j\ + 2,... . (2.26) 
Finally, defining 

Wfj = ex P (=Fpi) \jm) , (2.27) 
1 Mathematically oriented reader would prefer 

(A' A ® ls)|nA,ns) = nA\nA,riB) , (1a <8> B)|nA,ras) = n_B|nA,ns) , 
where 1a (Is) denotes the identity operator in "A-sector" ( "B-sector" ) . 
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one obtains 

mfj = Efj^fj , (2.28) 

with 

Ef m = 2?iwj ± ihrf(2m + 1) . (2.29) 

Let us emphasize that the eigenvectors corresponding to energies (2.29) cannot be 
normalized and should be considered as generalized eigenvectors not belonging to the 
Hilbert space of the problem. 



3 Canonical quantization of non-Hamiltonian systems 

As is well known any dynamical system may be regarded as a part of a larger Hamilto- 
nian system. Bateman's approach is based on adding to the primary system a time re- 
versed (dual) copy. Together they define a Hamiltonian system. There exists, however, 
a general approach to canonical quantization of non-Hamiltonian systems based on an 
old observation of Pontriagin [31]. Suppose we are given an arbitrary non-Hamiltonian 
system described by 

x = X(x) , (3.1) 

where X is a vector field on some configuration space Q. For simplicity assume that 
Q C R. N , that is, the system has N degrees of freedom. This system may be lifted to 
the Hamiltonian system on the phase space V = Q x R N as follows: one defines the 
Hamiltonian H : V — ► R by 

N 

ff(x, p) = p • X(x) = PlXlW , (3-2) 
l=i 

where (x, p) = (x±, . . . ,xn,Pi, ■ ■ ■ ,Pn) denote canonical coordinates on V . The corre- 
sponding Hamilton equations read as follows: 

x k = {x k ,H} = X k (x) , (3.3) 

Pk = {p k ,H} = -J2 Pl ?^-, (3.4) 

l=i Xk 

for k = 1, . . . , N. In the above formulae { , } denotes the standard Poisson bracket on 

WG}= £(«l*l (3 . 5) 

f^ydxkdpk dx k dp k J 

Clearly, the formulae (3.3) reproduce our initial dynamical system (3.1) on Q. The 
canonical quantization is now straightforward. Assuming the symmetric Weyl ordering 
one obtains the following formula for the quantum Hamiltonian 

^quantum = W ^ p,X,(x)^ , (3.6) 
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where W(/) denotes the Wigner-Weyl transform of a space-phase function / = /(x, p). 
Recall, that the Wigner-Weyl transform of / is defined as follows 



/ = W(/) = y da J dr/(<r,r)exp |i^(cj fc % + T fe p fc )j 



(3.7) 



where /(tr,r) denotes the Fourier transform of /(x, p). Clearly, -ff qU antum defines a 
Schrodinger system in L 2 (R N , c/x). 

Consider now a damped harmonic oscillator described by 

x + 27X + kx = . 

The above 2nd order equation may be rewritten as a dynamical system on M 2 

±1 = -jxi+ujx 2 , (3.8) 
X2 = —7x2 — ojxi , (3.9) 

with lo defined in (2.6). Clearly this system is not Hamiltonian if 7 / 0. However, 
applying the above Pontriagin procedure one arrives at the Hamiltonian system on M 4 
defined by the following damped harmonic oscillator Hamiltonian: 

H(x, p) = lo{pix 2 - P2X1) - 7(pixi + P2X2) ■ (3.10) 

The corresponding Hamilton equations of motion read 

x = Fx , p = -F T p , (3.11) 

where 

*= (::.;)• < 3 -> 

and F denotes the transposition of F. One may ask what is the relation between Bate- 
man's Hamiltonian (2.5) and that obtained via Pontriagin procedure (3.10). Surpris- 
ingly they are related by the following simple canonical transformation (x,y,p x ,p y ) — ► 
(xi,x 2 ,Pi,p 2 ): 

xi = , pi = -Vojy (3.13) 



x 2 = -\/ujx , p 2 = -~j= ■ (3.14) 

Assuming the symmetric Weyl ordering one obtains the following representation of the 
quantum Bateman's Hamiltonian (2.7) with 

H = uj(pix 2 -V2X1) , (3.15) 

and 

7 

Hi = -- (pixi + xipi +P2X2 + X2P2) ■ (3.16) 
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4 Spectral properties of the Hamiltonian 

4.1 Polar representation 

The formula (3.10) for H considerably simplifies in polar coordinates: 

x\ + 1x2 = re tLp . 
Denning the corresponding conjugate momenta 

xp 

P<p = L3 , p r = — , (4.1) 
r 

with L3 denoting 3rd component of L = x x p in M 3 , one finds 

H = -ujpp - jrp r . (4.2) 
The Hamilton equations in polar representation have the following simple form: 

<p = —uj , p v = , (4.3) 

and 

r = —jr , p r = 7p r . (4.4) 

The polar representation nicely shows that the Hamiltonian dynamics consists in pure 
oscillation in 93-sector and dissipation (pumping) in r-sector (p-sector) . In our opinion 
it is the most convenient representation to deal with . 
The quantization of (4.2) leads to (2.7) with 

d 

H = -ujp^p = iujh — , (4.5) 

and 

Hi = ijh (r-^- + lj = -7 (rp r - y^j , (4.6) 
where the radial momentum p r is defined by 

* = -"(I + 5f)- < 47) 

One easily finds the polar representation of the su(l, 1) generators: 
$ hf d 2 1 d 2 \ 1 2 

Jl = -AW 2 - + V 2 w)~Tn r ' (48) 



d 



■h = - 2 {ry r+ l), (4.9) 



together with the Casimir operator 



Note, that unitary evolution generated by H is given by 

(j(t) = e -iHt/h = e -iH t/h e -iHjt/h = e7 t exp (ut-^j exp (jt r-j^j , (4.12) 

and hence 

{U(t)ip)(r, if) = e 7 V(e 7 V, <p + ut) . (4.13) 
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4.2 Complete set of eigenvectors 

It is evident that H defines an unbounded operator in 7i = L 2 (R 2 , dxidx-z). It has 
continuous spectrum a(H) = (—00,00). To find the corresponding generalized eigen- 
vectors let us note that in polar representation the Hilbert space H of square integrable 
functions in M 2 factorizes as follows: 

L 2 (R 2 ,dx 1 dx 2 ) = L 2 ([0,2Tr),dip) ®L 2 (IR+, rdr) . (4.14) 

Therefore, the spectral problem splits into two separate problems in L 2 ([0, 2tt), dip) and 
L 2 (R + ,rdr). One easily finds 

HV lx = E lx V lx , (4.15) 

with 

Eix = h(lu + A7) . (4.16) 
The corresponding eigenvectors are defined by 

* lx {r,<p) =$ l (i P )R x (r) , (4.17) 

where 

$,(¥>):=-=, i = 0,±l,±2,... , (4.18) 

V Z7T 



and 

r -(a+i) 

R x (r) = —=- , A G R . (4.19) 

V Z7T 



Note, that G L 2 ([0, 27r), <i(^) whereas i?A does not belong to L 2 (R + ,rdr). 
One easily shows that the family ^i X satisfies 



and 



rZTT poo 

/ / Wx(r,p)*i>y(r,<p)rdrd<p = 5 lll 5(\-\'), (4.20) 
Jo Jo 

Y / Wx(r, <p) *ix(r', <p') dX = - 8(r - r')% - <p') . (4.21) 



l=—oo 

They imply the following resolution of identity 



_ / LA., 

1= £ / dA|*, A )<*| A |, (4.22) 



l=—oo ' 

and the spectral resolution of Hamiltonian 



00 /*0O 

H= Y, d\E lx \V lx ){V lx \, (4.23) 

i_ J— 00 



l=—oo 
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4.3 Feynman propagator 

Let us calculate the corresponding Feynman propagator 

K(x,t|x',t') = (x\U(t-t')\x.'} , (4.24) 
where U{t) = exp(—iHT/h). Using polar representation one finds 

00 POO 

K(r,tp,t\r', i p > ,t')= / e-^ T / h ^ lx {r^)Wx{r' \<p')d\ , (4.25) 

1 J —00 

with t = t — t' . Now, using (4.17) one obtains 

K(r,<p,t\r',<p',lf) = K 1 (r,t\r',t')K 2 (<p,t\y>',t') , (4.26) 
where the radial and azimuthal propagators are given by 

/CO 
e-^ T R x (r)R- x (r')dX, (4.27) 
-00 

and 

00 

K 2 (<p,t\<p',t') = Y, e-^MfWM, (4-28) 
{=— 00 

respectively. Finally, formulae (4.18) and (4.19) imply 

K 2 (<p, t\<p',t') = 6(<p' -<p-wr), (4.29) 

and 



1 1 f°° 

KAr.tW'J) = ——- / e iA(lnr'-lnr-7r) dX 



1 S(r' — re< T \ 
= _L ^(in r'-lnr- 7 r) = e 7r ^ -— -i . (4.30) 

Therefore, the time evolution is given by 

Y>t(r,¥>) = / / K{r,<p,ty,<ffJ = 0)Mr',<P f )r'dr'dtp' 
Jo Jo 

= e 7t Vo(e 7 V,^ + fjt) , (4.31) 
which perfectly agrees with (4.13). 

5 Analyticity and complex eigenvalues 

Now we are going to relate the energy eigenvectors *& n \ corresponding to the real 
spectrum E nX with the family of discrete complex eigenvalues of the Bateman's Hamil- 
tonian. Let us consider the distribution with A G C, i.e. for any test function 
<f>(r, if) 

tf, A (# = {^a} = / r- a Mr)dr, (5.1) 
J 
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where 

<M0 = ^ J Q e^<f>(r,p)dip. (5.2) 

Now, the analytical properties of depend upon the behavior of 4>i{r) at r = 0. A 
distribution r a acting on the space of smooth functions £(R+) 

POO 

S(R+) 3 f — / r a J(r) dr , (5.3) 
J o 

is well defined for all a £ C except the discrete family of points where it may have 
simple poles (see e.g. [32]). The location of poles depends upon the behavior of a test 
function / at r = 0. Assuming the most general expansion of f(r) 

f(r)=k + hr + f 2 r 2 + ... , (5.4) 

the poles are located at a = —1, —2, —3, .... However, 4>i(r) defined in (5.2) is much 
more regular. It can be observed (see Appendix B.) that 4>i{r) may be expanded at 
r = as follows: 

Ur) = ai r"! + a l+2 r^ 2 + a l+4 rl'l+ 4 + . . . . (5.5) 
Therefore, the poles that remain are located at 

Kl = -i(|Z| + 2n+l) , n = 0,1,2,... . (5.6) 
Moreover, the corresponding residues of are given by 

1 fnl 



Res^ 



where 



On the other hand 



A=A ni y/(\l\+2n)[ V2n 

j(_l)|J|+2n S (\l\+2n)^ 



(5.7) 



,nlK ^' V(KI+2n)! r y 



/A 



= V(|Z|+2n)!-fi=, (5.9) 

A=A„; V^TT 



where 



r |/|+2n 

^^TPW (5 ' 10> 



Now, the crucial observation is that satisfy 



and 



Jb |f£> = £ |f£> , (5.H) 



J 2 |f±)=±A(|Z| + 2n+l)|f±), (5.12) 



which proves that they define eigenvectors of H 

= Bftfc) , (5.13) 
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corresponding to complex eigenvalues 

E±=hwl±ih~j{\l\ + 2n + l). (5.14) 

The above formula is equivalent to the Bateman's spectrum (2.29) after the following 
identification 

3 = - r (5-15) 

and 

m =I(|/| + 2n) = |i|+n, (5.16) 
which reproduces condition (2.26). In terms of (jia^b) one has 

n A = I(|/| + /)+ n , (5.17) 

n B = \{\l\-l)+n. (5.18) 

We have therefore the following relation between \4>f m ) and |f^): 

l^> = l$,m- W >. ( 5 - 19 ) 
that is, |f±) defined in (5.8) and (5.10) may be regarded as a particular representation 

°fivf m >. 

Let us introduce two important classes of functions [33]: consider the space of 
complex functions / : C — ► C. A smooth function / = /(A) is in the Hardy class from 
above 7i\ (from below Tt 2 _) if /(A) is a boundary value of an analytic function in the 
upper, i.e. ImA > (lower, i.e. ImA < 0) half complex A-plane vanishing faster than 
any power of A at the upper (lower) semi-circle |A| — > oo. Now, define 

S_ = {0es\(V lx \<p) eH 2 _ } , (5.20) 

that is, (f> € 5- iff the complex function 

C9A — <^ A |0> eC, 

is in the Hardy class from below Ti 2 _. Equipped with this mathematical notion let us 
consider an arbitrary test function </> G S-. The resolution of identity (4.22) implies 

°° />0O 

<t>{r,f)= E / d\V lx (r,<p)(y lx \4>). (5.21) 

i J -co 

Now, since {*f?i\\4>) £ H-, we may close the integration contour along the lower semi- 
circle |A| — > oo (see Figure 1). 

Hence, due to the residue theorem one obtains 

oo oo 

<l>(r,<p) = -2m ^Res* a (r,^) <*ja|0> • (5-22) 



l=—oo n=0 



n 




Figure 1: Integration contour along the lower semi-circle for 1 = 2. 
Finally, using (5.7) and (5.9) one gets 

oo oo 
l=— oo n=0 



(5.23) 



We have proved, therefore, that the subspace 5_ C S C TL gives rise to the following 
resolution of identity 



1_ = 1 



E Eiua + 



(5.24) 



£=— oo n=0 

The same arguments lead us to the following spectral resolution of H restricted to 5_ : 



H- = H 



E E E nl \fnl)(fnl\ ' 



(5.25) 



l=—oo n=0 

with E 1 ^ defined in (5.14). Introducing the following family of operators 

^nl = \fnl )(fnl\ ' 

the spectral decompositions (5.24) and (5.25) may be rewritten as follows 



1L 



oo oo 

EE* 1 . 



nl i 



and 



Note, that 



l=— oo n=0 



oo oo 



H ~ = E E ^iZ • 



/=— oo n=0 



nl ' 



(5.26) 



(5.27) 



(5.28) 



(5.29) 



that is, the family P nl seems to play the role of the family of orthogonal projectors. 
Note, however, that are not hermitian. 
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6 Time reversal 

It was shown in [20] that Bateman's Hamiltonian is time reversal invariant 

T^HT = H , (6.1) 

where T denote the anti-unitary time reversal operator. Moreover, it turns out [20] 
that both Jo and J 2 satisfy 

tU t = Jo , rt j 2 r = j 2 . (6.2) 

Let us define 

Sia = Ttf , A . (6.3) 
In analogy with (4.22) and (4.23) one has the following resolution of identity 

00 POO 

1= Y, / dA|S IA )<S IA |, (6.4) 

i=-oo ,/ - 00 

and spectral resolution of the Hamiltonian 



^ rOG 

/ d\E lx \E lx )(E lx \. (6.5) 

J— 00 



H 

1^00 



Now, let us introduce another subspace S + in the space of test functions 

S + = [<pes\(E lx \ ( f>) en 2 + ], (6.6) 
that is, cj) £ 5+ iff the complex function 

C 9 A — > (E^) eC, 
is in the Hardy class from above H\. It is easy to show that 

S + = T(<S_), (6.7) 

and vice versa 

5_ = T(S + ) . (6.8) 
Indeed, if G S- then ( ^i X \4>) £ Ti^_. One has therefore 

( E lx \T<t> ) = < 0|Tt E lx ) = g , (6.9) 

which implies that T(f> 6 5+. 2 Moreover 

5_ n S+ = {0} . (6.10) 

2 In the above formulae we have used 
which holds for any anti-linear operator A. 
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To prove this property let us assume that <p G 5_ fl 5 + . Since (f> G S + , one has 
(H JA |0) G H+. However 



(~ lx \ct>) = (cf>\T* lx ) = (y lx \T*cf>} eH 



+ 



(6.11) 



On the other hand 7"t G S_ and hence ( ) G Therefore, (*; A |Tt</>) G 

H 2 . fl which means that (vE^IT </>) is en entire function vanishing on the circle 
|A| — ► oo. However, any entire function is necessarily bounded and hence such <f> 
which belongs both to 5_ and S + does not exist. 

Now, take any test function <p G S+. Formula (6.4) implies 



°° /"OO 

<K r >¥0 = E / dASiA(r,¥>)<Ej A |0). 

J_ J -co 



(6.12) 



Let us continue the eigenvectors "Ei\ for the complex A plane. They have simple poles 
at A = —X n i with A„; defined in (5.6). The corresponding residues of follows from 
(5.7) 

1 r r nl 



Res S 



iA 



Moreover, 



-/A 



a=-a„ ; ^(j/j +2n)! v 7 ^ ' 
= V(|/|+2n)!^fi. 

A=— A n ; V^VT 



(6.13) 



(6.14) 



Now, since (E n \\cf)) G H+, we may close the integration contour in (6.12) along the 
upper semi-circle |A| — > 00. The residue theorem implies 



00 00 



(r,(p)=2m E Res 



l=— 00 n=0 

Finally, using (6.13) and (6.14) one gets 



A — — A^ 



(Z lx \cf>) 



A— — A^ 



(6.15) 



Z=— oo n=0 

and hence it implies the following resolution of identity on S + : 



1 



5+ 



E E r iuai^ f 



)7l=0 



Now, the formula (5.12) together with (6.2) gives 



■fcr|f±> = =F-(|l| + 2n+l)T|ft> 



(6.16) 



(6.17) 



(6.18) 



and hence one deduces the following relations between |f^) and time reversed T\f^ 



n + nl ) 



iq , 



o ta nl 



(6.19) 
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where a n \ are arbitrary (n, Z)-depended phases. It should be stressed that these phases 
are physically irrelevant. Actually, one may redefine |f^) in (5.8) and (5.10) such that 
these additional phase factors disappear from (6.19). Let us observe that 

T 2 \ti) = \ti) , (6-20) 
irrespective of a n [. Taking into account (6.19) one obtains from (6.17) 

oo oo 
+ Z=-oon=0 

The same arguments lead us to the following spectral resolution of H 

oo oo 

A + = A \ S = E E^oy, (6.22) 

l=— oo n=0 

with E^ t defined in (5.14). Finally, introducing 

P:i = Kl)(rj = {PJ , (6.23) 

with P~ t defined in (5.26), the spectral decompositions (6.21) and (6.22) may be rewrit- 
ten as follows 

oo oo 
/=— oo n=0 

and 

oo oo 
Z=—oo n=0 

7 Resonances and dissipation 

What is the physical meaning of the complex eigenvalues E nl ? To answer this question 
let us consider the resolvent operator of the Bateman's Hamiltonian 

R{H,z) = (H- z)- 1 . (7.1) 

Using the family of eigenfunctions ) one has 

00 r°° d\ 

R(H,z) = Y, / p 1 I*iaX*ia|, (7-2) 

with Ei\ defined in (4.16). Now, using the same technique as in Section 5 one easily 
finds 



oo oo 



r_(*) = roML = E EWtt^' ( 7 - 3 ) 

^ /=-oon=0 ^nl Z 

with P~ t defined in (5.26). This shows that E~ t constitute poles of the resolvent operator 
on S- . In the same way using the family | Ei\ ) 

R(H,z) = £ / -P— \E lx )(E lx \, (7.4) 
, — J-oo Ea - Z 



l=—oo 
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ijh(\l\ + 2n+ 1) 



hujl 



Figure 2: A closed curve 7~t on a complex energy plane. 



one finds 



R+(z) = R(H,z) 



s+ 



EE 



E + — z 

l=-oon=0 nl 



P. 



nl ' 



(7.5) 



which shows that Ej~, constitute poles of the resolvent operator on 5+. As is well 
known [29] the poles of the resolvent operator correspond to resonant states. Hence, 
the complex eigenvalues E^ may be interpreted as resonances of the Bateman's Hamil- 
tonian. Note that due to the Cauchy theorem operators may be represented by the 
following integrals 



1 

2ni 



"hil 



(7.6) 



where 7 , is any (clockwise) closed curve which encircles a single pole z = E , (see 
Figure 2). 

Finally, let us turn to the evolution generated by the Bateman's Hamiltonian. 
Clearly, 

R B t — ► U(t) = exp(-iHt/fi) , 

defines a group of unitary operators on the Hilbert space L 2 (M 2 ). Now, it is easy to 
see that if tp- G 5_, then U(t)ip- belongs to 5_ only if t > 0. Similarly, if ip + £ S + , 
then U(t)ip + belongs to 5+ only if t < 0. Therefore, we have two natural semigroups 



and 
where 

One has 



U-{t) : S 
U+(t) : 5. 



and 



for t > , 

for i < , 
#+(t) = £>(i) 



(7.7) 

(7.8) 
(7.9) 

(7.10) 



n=0 
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for t > 0, and 

oo oo 
l=— oo n=0 

for £ < 0. It should be clear that these two semigroups are related by the time reversal 
operator T: indeed formulae (6.19) imply 

TP'ai^ = Pni ^d TP+Tt = p-, (7.12) 

and hence 

oo oo 

l=—oo n=0 

oo oo 

l=— oo n=0 

for t > 0. Similarly, one finds 

rt)+(t)rt = ^_(-t) , (7.i4) 

for t < 0. We have shown that perfectly reversible quantum dynamics f7 (t) on the 
full Hilbert space L 2 (K. 2 ) is no longer reversible when restricted to the subspaces 5_ 
and S + . This effective irreversibility is caused by the presence of resonant states | f ^ ) 
corresponding to complex eigenvalues E^. 

8 Conclusions 

In this paper we have studied the spectral properties of the Bateman Hamiltonian. It 
was shown that the complex eigenvalues E^ m given by (2.29) corresponds to the poles of 
the resolvent operator R(H,z) = (H — z)" 1 . Therefore, the corresponding generalized 
eigenvectors may be interpreted as resonant states of the Bateman dual system. It 
proves that dissipation and irreversibility is caused by the presence of resonances. 

From the mathematical point of view the Bateman system gives rise to the so called 
Gel'fand triplet or rigged Hilbert space [26, 27] (see also [28, 34]). A Gel'fand triplet 
(rigged Hilbert space) is a collection of spaces 

$CHC$', (8.1) 

where H is a Hilbert space, $ its dense subspace and <&' is the dual space of continuous 
linear functionals on <!>. Note, that elements from do not belong to TC. This is 
a typical situation when one deals with the continuum spectrum. The corresponding 
generalized eigenvectors are no longer elements from the system Hilbert space. They 
are elements from the dual space i.e. distributions acting on <I> [35, 36]. In our case 
we have two natural Gel'fand triplets: 

S- C L 2 (R 2 ) C S'_ , (8.2) 



17 



and 



S + C L 2 (M 2 ) C 5; 



(8.3) 



The first triplet corresponds to the forward dynamics U- and the second one corre- 
sponds to the backward semigroup U + . A similar analysis based on rigged Hilbert 
space approach was performed in [24, 25] for a toy model damped system defned by 
x = — 7X. 



Appendix A 

Let us briefly sketch calculations leading to (5.7) and (5.9). We introduce a distribution 
acting on a test function 0(r, ip) as an antilinear functional defined by the integral 



tf, A ($ = (<f>\Vix) 



1 

27 



e -a Vr -i\-i- 



{r,<p)dS= I r _i Vi(r)dr, (A.4) 
/o 



where X £ C , dS = rdrdip, and 4>i{r) is given by (5.2). Expanding (j>i(r) in the power 
series and rewriting the last integral as 



f 

Jo 



r lX (j)i(r)dr 



+ 
+ 



/ 

Jo 

roc 
Jl 

f 

Jo 



-iX 



^(r)-^(0)-r^(0) 



r l 1 -(i-i) 



(0) 



-iX~T 



r)dr 



dr 

(A.5) 



,-iX 



J (O)+r^(O) + ... + 77 ^T^ >) 



(7^1)! w 



one can observe that the first two summands are regular for all A G C. The last integral, 
however, equals to 



■'0 I — n 



(0) 



1 



k=0 



k=0 



k\ k-iX+1 



(A.6) 



and has simple poles in A& = —i(k + 1), k = 0, 1, . . . , I — 1. Moreover, one can read 
from (A.6) that 

A=-i(fe+i) A;! 



Res (<6|* 



/A; 



Finally, using 



and 



4 k \o) = ^j\- il ^ {k \o,<p)d v 



<f> {k \o, ( p) = (-l) k r S -^M (f) (r,<p)rdr 
Jo r 



we get 



0z (0) - 



2vr 7r2 



/ 

iff 



r 



(A.7) 
(A.8) 
(A.9) 
(A.10) 
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But due to (5.5) in the case investigated here k = \l\ + 2n, hence the poles are located 
at X n i = + 2n + 1) and 



Res 



nil 



A=A ni y/(\l\+2n)\ V 7 ^ ' 



(A.ll) 



where f nl is a distribution given by (5.8). 

The conjugate distribution is defined as 



V lx (<j>) = (cf)\y lx ) 



1 



(r, 



r* A ^;(r)dr 



2vr 7 R 2 j 
and it is regular in A = A n j. The poles of (4>\^ix) are located at A = A n ^. Hence 

ifS 



(A.12) 



/A/ 



A — 



V(l'l+2n)! 



(A.13) 



where is a distribution given by (5.10). 

Appendix B 

Let us briefly proof that 4>i(r) given by (5.2) has a power series expansion (5.5) starting 
from rl'L Supposing that (f>(xi,X2) is an analytic function 



ki,k 2 UX 1 aX 2 



.T. 



k 2 



(B.l) 



in cartesian coordinates (xi, X2) it is obvious that in polar (r, (^-coordinates one obtains 
the following expansion for <j>i{r): 

1 f 2lT 

<M r ) = — / e~ d(p 4>(r cos (p,r sin ip) dip 
2vr Jo 

= 7T e"^ 0(0,0)+ V A fel , fe2 r fel+fe2 (cos^) fel (sm^) fe2 ^ 

27T Jo L trL J 



fci,fc 2 



where 



-!- J] ^ llfc2 r fcl+fc2 / e-^(cos^) fcl (sm^) fe2 ^, 
d kl+k2 (p(0,0) 



(B.2) 



A 



dx kl dx k 2 2 



stand for derivatives of 0(xi,X2) in (0,0). Now, the question is: for which values of 
k = k\ + /c2 the sum in (B.2) does not vanish? Clearly, it should be 



/•27T 

L ' 



e~^(cos ip) kl (sin iff 2 dip^O. 



\ki 



k 2 



(B.3) 
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Using the Newton expansions 

mi=0 V m l/ 

1 &2 & 1 ^2 ^ f Jc \ 

(sinv9) fc2 = (-) '(e^-e"^) * = (— ) ' ^ 2 J {-\) k ^-^ e im ^ e -i{k 2 -m 2 ) v ^ 

1 m 2 =0 V m 2/ 

one can rewrite (B.3) as 
hence (B.4) will not vanish iff 

l + k — 2{m\ + m 2 ) = . (B.5) 

Clearly, 

< mi < fci , < m 2 < k 2 . (B.6) 
Now, let I < 0, so Z = — \l\ and 

fc = |/| + 2(mi + m 2 ) = |/| + 2n , (B.7) 

where n = mi + m 2 > 0. Due to (B.6), in order to satisfy (B.7) it should be k > \l\. 
On the other hand, if / > 0, then 

k = -l + 2(m 1 +m 2 )<-l + 2k, (B.8) 

because of (B.6) and finally k > I = \l\. Note that in this case k = I + 2(m\ + m 2 — I), 
where m\ + m 2 — I = n > 0. 

As a result we obtained that for a given n the lowest power of r in expansion (B.2) 
is k\ + A; 2 = Moreover, in both cases 

k 1 + k 2 = k = \l\+2n, n = 0,1,2,... (B.9) 

holds. 
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